Dynamical turbulent flow on the Galton board with friction 
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We study numerically and analytically the dynamics of particles on the Galton board, a regular 
lattice of disc scatters, in the presence of a constant external force and friction. It is shown that 
under certain conditions friction leads to the appearance of a strange chaotic attractor in an initially 
conservative Hamiltonian system. In this regime the particle flow becomes turbulent and its average 
velocity depends in nontrivial manner on friction and other system parameters. We discuss the 
applications of these results to the transport properties of suspended particles in a laminar viscous 
flow streaming through scatters. 
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It is well known that dissipation can lead to the appear- 
ance of strange chaotic attractors in nonlinear nonau- 
tonomous dynamical systems In this case the en- 

ergy dissipation is compensated by an external energy 
flow so that stationary chaotic oscillations set in on the 
attractor. Such an energy flow is absent in the Hamilto- 
nian conservative systems and therefore the introduction 
of dissipation or friction is expected to drive the system to 
simple fixed points in the phase space. This rather gen- 
eral expectation is surly true if the system phase space 
is bounded. However a much richer situation appears 
in the case of unbounded space, where unexpectedly a 
strange attractor can be induced by dissipation in an 
originally conservative system. To investigate this situ- 
ation we study the dynamics of particles on the Galton 
board in the presence of a constant external field and 
friction. This board, introduced by Galton in 1889 ||, 
represents a triangular lattice of rigid discs with which 
particles collide elastically. For the case of free particle 
motion, the collisions with discs make the dynamics com- 
pletely chaotic on the energy surface as it was shown by 
Sinai [Q. In this paper we study how the dynamics is 
changed if the particles are affected by a friction force 
Ff which is assumed to be proportional to the particle 
velocity v and directed against it: Ff — —7V. In the free 
space without discs, an external field with an effective 
force f creates a stationary particle flow with the veloc- 
ity Vf = f/7. All perturbations decay to this flow with 
the rate proportional to 7 so that this laminar flow can be 
considered as a simple attractor. It is interesting to un- 
derstand how this flow will be affected by the scattering 
on regularly placed discs. As we will show latter on, in 
this case the flow can become turbulent and chaotic, and 
demonstrate many unusual properties (see Fig.l). The 
Hamiltonian dynamics of the model has been studied in 
great detail (see || and Refs. therein) however, as for 
our knowledge, the effects of friction were not addressed 
up to now. The investigation of the dynamics in this 
regime is interesting from the viewpoint of general phys- 



ical grounds. But we also note that physically such a 
situation can correspond to the motion of small particles 
in a viscous liquid flow which laminarly streams through 
large scatters that in our case have the form of discs. The 
same model can describe the electron dynamics in anti- 
dot superlattice which has been experimentally realized 
in semiconductor heterostructures [|| . In such structures 
the effects of classical chaos play an important role 
and the effects of friction we discuss here can appear for 
relatively strong electric fields. 




FIG. 1. Chaotic flow on the Galton board. Here the dis- 
tance between discs is R — 2.24, f = (-0.5, -0.5) and 7 = 0.1. 
Initially 100 particles are distributed homogeneously along a 
straight line segment in the upper right corner, their color 
homogeneously changes from red to green along this segment. 

To study the dynamics of this model, we fix the disc 
radius a = 1, and the mass of the particles m — 1 so 
that the system is characterized only by three parame- 
ters: the distance R between discs packed into equilateral 
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triangles all over the (x, y) plane; the friction coefficient 7 
and the external force of strength /. For R< R c = 4/V3 
there is no straight line, by which a particle can cross the 
whole plane without collisions (at / = 0, 7 = 0), and we 
start the discussion from this case. The particle dynam- 
ics is simulated numerically by using the exact solution of 
Newton equations between collisions, and by determin- 
ing the collision points with the Newton algorithm. This 
way the trajectories are computed with high precision, 
and for example for 7 = the total energy is conserved 
with a relative precision of 1CT 14 . A typical example of 
the chaotic flow formed by an ensemble of particle trajec- 
tories is shown in Fig. 1. To illustrate the mixing prop- 
erties of the flow we attributed a color to each trajectory 
that allows to follow their interpenetration and spread- 
ing. This figure shows that there is a certain penetration 
depth of one color into another, however this depth is 
finite since on average the initial color repartition is still 
visible. 




-9700 1 1 1 1 

-10000 -9990 X 9980 

FIG. 2. Example of a single trajectory for the case of 
Fig. 1 shown on a small (main figure) and large (inset) scales. 
The drift velocity of the flow is Vf w 0.13. 

The properties of color penetration can be understood 
from the analysis of single trajectory dynamics. Such a 
typical example is presented in Fig. 2. It shows that tra- 
jectory flows with an average constant velocity Vf in the 
direction of the external force. This drift velocity is con- 
stant only on average since on a smaller scale the particle 
moves chaotically between scatters (see Fig. 2). Even if 
the average flow is regular, on a small scale the particle 
dynamics is chaotic and since the dynamics is dissipa- 
tive this means that the motion develops on a strange 
chaotic attractor. This attractor is induced by dissipa- 
tion in a conservative Hamiltonian system. For the case 
of Fig. 2 the drift velocity is relatively large and the par- 
ticle does not have enough time to move around many 
scatters in the direction perpendicular to the flow. As 



a result the penetration depth for color mixing is not 
very large. However for smaller friction, the drift veloc- 
ity becomes smaller and the penetration depth increases 
so that the particle makes many turns around discs as 
it is shown in Fig. 3. Surprisingly, the drift velocity v/ 
increases with the growth of the friction coefficient 7. 
This dependence is opposite to the case without scatters 
where Vf drops with the 7 growth. 
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FIG. 3. Single trajectory for the same parameters as in 
Figs. 1,2 but with smaller friction 7 = 0.004; on a large scale 
the particle moves with a drift velocity Vf « 0.05 in the di- 
rection of external force. 
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FIG. 4. Dependence of the flow velocity Vf on the friction 
coefficient 7 and external field strength / for R = 2.24, f/|f| 
= (-1, -1) and 0.001 < 7 < 0.4; 0.5 < / < 32. Numerical 
results are shown by circles, the straight line shows the de- 
pendence Vf oc (7/) 1 ^ 3 - Inset shows the dependence of Vf on 
AR = R - 2 for 7 = 0.1 and f = (-0.5, -0.5): points give 
numerical data, lines show the slopes 0.5 and 1. 

The dependence of the drift velocity Vf on 7 and on the 
external force / = |f | is obtained by averaging over many 
trajectories. Each trajectory is followed for a sufficiently 
long time so that the distance from the origin becomes 
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larger than a few thousands units. Up to fluctuations, the 
flow is always directed along the applied force f . The re- 
sulting dependence is shown in Fig. 4 according to which 
Vf oc (7/) 1 / 3 . To understand the origin of this unex- 
pected dependence on parameters we also compute the 
average velocity square v 2 =< v 2 + vr. > and determine 
how it varies with / and 7. The data presented in Fig. 5 
give a dependence of the form: v 2 oc f 4 ^ 3 /^ 2 ^ 3 - We stress 
that Vf and v 2 are independent of initial velocity of the 
particle. 

To explain the observed scaling relations for the dy- 
namical turbulent flow on the Galton board we give the 
following theoretical arguments. In the regime with weak 
friction the particles start to diffuse among discs in a 
chaotic manner with the diffusion rate D = vl/2 where v 
is the velocity of the particle and I is the mean free path 
H . For R ~ 1 we have I ~ R ~ 1 while the dependence 
of I on R will be discussed in more detail latter. During 
the dissipative time scale r 7 = m/7 this diffusion leads to 
the particle displacement Ar ~ yj Dm/j along the direc- 
tion of applied field. Here, for generality we introduced 
the particle mass m. This gives the change in the po- 
tential energy U ~ /Ar ~ /v/flm/7. In the stationary 
regime at time t 3> t 7 , this potential energy should be 
comparable with the kinetic energy of the particle so that 
U <~ mv 2 . As a result we find that the average velocity 
square is 

v 2 ~ fl\lhmfl 3 , (1) 
that is in agreement with the data of Fig. 5. 




FIG. 5. Dependence of average v 2 on / and 7 

for R = 2.24, f/|f| = (-1, -1) and parameter range 
0.001 < 7 < 0.4; 0.5 < / < 32. Circles show numerical data, 
the straight line gives the slope 1. 

The relation (Q) allows to determine the mobility fx 
which relates the velocity of the flow with the applied 
force: Vf = /1 f. Indeed during the time r c between 
collisions the particle is accelerated by force that gives 
the average drift velocity Vf — fr c /m and fx — r c /m. 



Since the dynamics is chaotic the direction of velocity is 
changed randomly after each collision so that vt is ac- 
cumulated only between collisions. The time t c is deter- 
mined by the mean free path I and the average velocity 
v : Tc =l/v = 2D/v 2 and ,u = D/(mv 2 /2). In fact 
this is the Einstein relation according to which the mo- 
bility is given by ratio of the diffusion rate to the average 
kinetic energy (temperature) 0. Together with ([!]) this 
gives the drift velocity of the chaotic flow: 

Vf-l^Wm*) 1 '* , (2) 

which is in a good agreement with the numerical results 
in Fig. 4. On average the particle moves with the velocity 
Vf in the direction of the force (see an example in the 
inset of Fig. 2). The amplitude of fluctuations around 
this direction is Ar ~ y/Dm/j which also determines 
the color mixing depth (see Fig. 1). 

The relations ([!]) , (||) allow to estimate the value of the 
Lyapunov exponent A. Indeed, the particle moves with a 
typical velocity v and as in the case of the Sinai billiard 
with I ~ R ~ 1 we have A ~ v/l. Therefore in the regime 
when: 

7 < 7c ~ V / 7™A , (3) 

the value of A is much larger than the dissipation rate 
7/m. As a result for 7 <C 7 C the strange attractor is fat 
and its fractal dimension is close to the maximal dimen- 
sion 4, which is determined by the number of degrees of 
freedom (we remind that contrary to the nondissipative 
case the energy is not conserved). 

The physical meaning of the critical "f c border is the 
following: for 7 » 7 C the dissipation time r 7 becomes 
much shorter than the time between collisions r c . In this 
case the dissipation dominates chaos and the strange at- 
tractor degenerates into a simple attractor. Among such 
simple attractors we observed trajectories that move reg- 
ularly under some angle to the applied force but more 
often trajectories that simply drop to a fixed point on 
a disc where the applied force is perpendicular to the 
disc border. Of course the numerical simulations can not 
guarantee that the chaotic attractor at 7 < 7 C will not 
degenerate to a simple one at very long times. However 
all our simulations performed with high computer preci- 
sion give no indications in this direction even when we 
followed trajectories up to very far distances from the 
origin (bigger than 10 5 length units). 

The equation (Q), which determines the border for 
disappearance of chaotic flow, is written for the case 
I ~ R ~ 1. Therefore, it should be modified when R 3> 1 
and AR = R - 2 < 1. Indeed for R > 1 the mean 
free path I is given by / = l/(<rn) where a = 2a = 2 is 
the disc cross section and n — 2/(y/3R 2 ) is the density 
of scatters. Then the time between collisions is l/v and 
from equation (^) we obtain j c ~ mv/l ~ \J fma/ R. In 
the other limit AR < 1 we may assume that D ~ vAR, 
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I ~ Ai? and the escape time from a cell between three 
discs is r esc ~ a 2 /D. Therefore the dissipation be- 
comes dominant and the chaotic attractor disappears at 
7 > 7 C - m/r esc ~ ^/mf/a 3 Ai?. 

The above changes in the mean free path I also af- 
fect the drift velocity of the flow through the relation 
(^). Indeed for Ai? « 1 we expect I ~ Ai? that gives 
vj oc Ai? 2 / 3 . This dependence is close to the numer- 
ical data shown in Fig. 4 (inset) even if the numerical 
value of the exponent is approximated better by 0.5. In 
the other limit i? ^> 1, we have I ~ i? 2 /a that gives 
Vf oc i? 4 / 3 . This power dependence is in a satisfactory 
agreement with our numerical data (inset in Fig. 4) al- 
though the numerical value of the exponent is closer to 
1. We attribute these small deviations in the exponent 
values to a restricted interval of variation in i?. Actually 
we can not use very large/small values of Ai? since in 
these limits the value of 7 becomes comparable with j c 
and chaotic attractor disappears. We note that according 
to our data (see Fig. 4) the strange attractor exists even 
in the case i? > R c = 4/\/3 when at / = 0, 7 = there 
are straight trajectories crossing the whole plane without 
collision. Apparently the contribution of these orbits is 
not significant if 7 > and f is not directed along these 
lines. 

Let us now discuss a few consequences of the ob- 
tained results. The introduction of scatters qualitatively 
changes the flow properties: instead of a laminar flow 
we obtain a turbulent/chaotic one. The flow rate Vf de- 
pends on the applied force / and 7 in a nonlinear way 
given by equation (^|). Contrary to the laminar case 
where Vf = //7 m presence of scatters the flow veloc- 
ity depends also on the mass of particles m. In fact light 
particles stream with higher velocity. This property of 
chaotic flow can be used for mass separation of particles 
suspended in viscous flows, for example isotopes. It is 
also interesting to note that in the derivation of equations 
(|l|), (H), (H) we didn't use any specific properties of the 
disc distribution in the plane. Therefore these results re- 
main also valid for a random distribution of disc scatters 
in the plane. We wonder if some rigorous mathematical 
results can be obtained in the limit of weak friction to 
prove the existence of a chaotic attractor for a regular or 
random distribution of discs. 

We think that the obtained results represent certain in- 
terest for transport properties of particles suspended in 
a viscous flow streaming through a system of large scat- 
ters. Indeed, when a liquid streams laminar ly through 
scatters with the velocity v s the flow of suspended par- 
ticles is characterized by the relations (0), (||), (||) with 
feff = v slejf where Jeff is an effective friction created 
by the viscosity of the liquid and f e f / is the resulting ef- 
fective force acting on the suspended particles. Then for 
weak friction the drift velocity of particles Vf is smaller 
than the velocity of the liquid v s . Such kind of transport 



can be studied experimentally with viscous liquids and 
its investigation can contribute to a better understand- 
ing of the interplay between dissipation, turbulence and 
chaos. 
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